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Abstract 

It has been well-known that for two-way contingency tables with fixed row sums 
and column sums the set of square- free moves of degree two forms a Markov basis. 
However when we impose an additional constraint that the sum of a subtable is 
also fixed, then these moves do not necessarily form a Markov basis. Thus, in this 
paper, we show a necessary and sufficient condition on a subtable so that the set of 
square-free moves of degree two forms a Markov basis. 

1 Introduction 



Since Sturmfels 11996 



1998 showed that a set of bino- 



and iDiaconis and Sturmfels 
mial generators of a toric ideal for a statistical model of discrete exponential fami- 
lies is equivalent to a Markov basis and initiated Markov chain Monte Carlo approach 
based on a Grobner basis computation for testing statistical fitting of the given model, 
many researchers have extensively studied the structure of Markov b a ses for model s 



in computational algebrai c statistics (e . g. iHosten and SullivantI 20021 ]. iDobral 2003 



Dobra and SuUivantI |2004| . lOeiger et all (2006|). 

In this article we consider Markov bases for two-way contingency tables with fixed 
row sums, column sums and an additional constraint that the sum of a subtable is also 
fixed. We call this problem a two-way subtable sum problem. From statistical viewpoint 



1 



this problem is motivated b y a block interaction model or a two-way change-point model 



proposed by iHirotsd | l997l| . which has been studied from both theoretical and practical 
viewpoint flNinomiyal 20041 ]) and has important applications to dose-response clinical trials 
with ordered categorical responses. 

Our model also relates to the quasi-independence model for incomplete two-way con- 



tingen cy tables which contain some structural zeros (lAoki and Takemural |2005| , iRapallo 
Essentially the same problem has b een studied in de t ail from algebraic view- 



Ohsugi and Hibi 



oint in a series of papers by Oh sugi and Hibi (lOhsugi and Hibil |1999a 
1999bt . [Qh"sugi and Hibil [200^ ). 

It has been well-known that for two-way contingency tables with fixed row sums and 
column sums the set of square-free moves of degree two forms a Markov basis. However 
when we impose an additional constraint that the sum of a subtable is also fixed, then 
these moves do not necessarily form a Markov basis. 

Example 1. Suppose we have a 3 x 3 table with the following cell counts. 



7 


5 


1 


5 


10 


6 


2 


6 


8 



// we fix the row sums (13,21, 16) and column sums (14,21, 15) , and also if we fix the 
sum of two cells at (1, 1) and (2, 1)(7 -|- 5 = 12 in this example), a Markov basis consists 
of square-free moves of degree two. However, if we fix the sum of two cells at (1,1) and 
(2,2)(7-|- 10 = 17 in this example), then a Markov basis contains non-square-free moves 
such as 



1 


1 


-2 


-1 


-1 


2 












In this paper we show a necessary and sufficient condition on a subtable so that a 
corresponding Markov basis consists of square-free moves of degree two. The results here 
may give some insights into Markov bases for statistical models for general multi-way 
tables with various patterns of statistical interaction effects. 

Because of the equivalence between a Markov basis and a set of binomial generators of 
a toric ideal, a theory of this paper can be entirely translated and developed in an algebraic 
framework. However, in this paper we make extensive use of pictorial representations of 
tables and moves . Therefore we prefer to develop our theory using tables and moves. See 



Aoki et al.l |2005l | for a discussion of this equivalence. 

This paper is organized as follows: In section [2l we describe our problem and summa- 
rize some preliminary facts. Section [3] gives a necessary and sufficient condition that a 
Markov basis consists of square-free moves of degree two. We end this paper with some 
concluding remarks in Section m 
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2 Preliminaries 



2.1 Subtable sum problem and its Markov bases 

Let N = {0, 1,2,.. .} and let X = {xy}, Xij e N, i = 1, . . . , R, j = 1, . . . ,C , he an R x C 
table with nonnegative integer entries. Let I = \ 1<J< C}. Using 

statistical terminology, we call X a contingency table and I the set of cells. 
Denote the row sums and column sums of X by 



C R 

Xi , 

=1 1=1 



^ ^ Xjj, i 1, . . . , R, x^j ^ ^ Xij, j 1, . . . , C. 



Let S' be a subset of X. Define the subtable sum x{S) by 



{i,i)65 



Xij. 



Denote the set of row sums, column sums and x{S) by 

b = . . . , xr+, . . . , x+c, x{S)}. 

For 5 = or S* = X, we have a;(0) = or x(I) = degX := — 
these cases x{S) is redundant and our problem reduces to a problem concerning tables 
with fixed row sums and column sums. Therefore in the rest of this paper, we consider 
5" which is a non-empty proper subset of X. Also note that x^S"^) = degX — x{S), where 
5''^ is the complement of 5". Therefore fixing x{S) is equivalent to fixing x{S'^). 

We consider b as a column vector with dimension R + C + 1. We also order the 
elements of X with respect to a lexicographic order and regard X as a column vector 
with dimension |X|. Then the relation between X and b is written by 

AsX = b. (1) 

Here As is an (i? + C + 1) x |X| matrix consisting of O's and I's. The set of columns 
of As is a configuration defining a toric ideal Ias- In this paper we simply call As the 
configuration for S. The set of tables X G N"^ satisfying ([1]) is called the fiber for 6 and is 
denoted by J^{b). 

An R X C integer array B = (ij)eJ satisfying 

AsB = (2) 

is called a move for the configuration As- Let 

Ms = {B\ AsB = 0} 

denote the set of moves for As- Let i3 C A^s be a subset of M.s- Note that if S is a 
move then —B is a mo v e. W e call B sign-invariant if i? G i3 ^ —B G B. According to 



Diaconis and Sturmfelsl |1998| . a Markov basis for As is equivalent to a set of binomial 



generators of the corresponding toric ideal for I^g and defined as follows. 
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Definition 1. A Markov basis for As is a sign-invariant finite set of moves B = {Bi, . . . , B^} 
C M.S such that, for any b and X , Y E J^{b), there exist a > 0, i?^^, . . . , Bt^ G B such 
that 

a a 

Y = X + and F = X + ^ Bj^, G J'(b) forl<a<a. 

s=l s=l 

In this paper, for simplifying notation and without loss of generality, we only consider 
sign-invariant sets of moves as Markov bases. 

For i ^ i' and j ^ j', consider the square- free move of degree two with +1 at cells 
ihj), i^'J') and -1 at cells (i,/) and (i'J) : 

3 j' 
i 1 -1 
i' -1 1 



For simplicity we call this a basic move and and denote it by 

B{t,t';j,f) = {t,j){t',f)-{t,f){t',j). 
It is well-known that the set of all basic moves 



Bo = J, /) I (z, j) G J, (z',/) el,t^ i\j ^ /} 

forms a unique minimal Markov basis for A%, i.e. the problem concerning tables with fixed 
rows sums and column sums. If B{i,i'; j, j') G A4s, we call it a basic move for S. Define 

Bo{S) = BonMs 

which is the set of all basic moves for S. Note that Bo{S) coincides with the set of 
square- free moves of degree two for As, since the row sums and columns sums are fixed. 

As clarified in Section[3l Bo{S) does not always form a Markov basis for As- In Section 
[31 we derive a necessary and sufficient condition on S that Bq{S) is a Markov basis. 



2.2 Reduction of Li-norm of a move and Markov bases 



In proving t hat Bn(S) is a Markov ba sis f or a given S, we employ th e norm-reduction 

Suppose that we 



argument of iTakemura and Aokil |2005l | and lAoki and Takemural [2003 
have two tables X and Y in the same fiber JF. Denote 



X -Y = {xij - yij}(ij)(.i 



and define the Li-norm of X — F by ||X — V||i = Yl 
\\X — Y\\i can be reduced (in several steps) by ^0(5") as follows. 



(i,i)6J \^i3 



We define that 
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Definition 2. For X in the same fiber T , we say that \\X — Y\\i can be reduced by 
Bo{S) if there exist r+ > 0, > 0, r+ + > 0, and sequences of moves Bf G Bq{S), 
t = 1, . . . , , and G Bq{S), t = 1, . . . ,t~ , satisfying 



\X -Y + Y,Bt + J2Br\U < \\X -Y\\ 

t=i t=i 

t' 

X + 5^i?+G^, for r' = l,...,r+, 

t=i 

t' 

Y-Y.BteJ', for r' = l,...,r-. 



1) 



(3) 



t=i 



In Takemura and Aoki |2005| we have mainly considered the case that ||X — y||i 
can be reduced in o ne st ep: r"*" + r~ = 1. However as discussed in Section 4.2 of 
Takemura and Aokil 2005l |. it is clear that Bq{S) is a Markov basis for As if for every 
fiber J^{b) and for every X ^ Y m. J^{b), \\X — Y\\i can always be reduced by Bq{S). 



Here the number of steps r"' 



r 



needed to reduce ||X — can depend on X and Y . 



Therefore we consider a condit i on th at \\X — can be reduced by Bq{S). 



As in lAoki and Takemural [2003l |. we look at the patterns of the signs of X 
Suppose that X — Y has the pattern of signs as in Figure [T]-(i)- This means 



Y. 



and the signs of Xij — yij and Xj/j/ — yiiji are arbitrary. Henceforth let * represent that the 
sign of the cell is arbitrary as in Figure [H Because Xi'j > 0, Xij> > 0, we have 

Vi'j > 0, yif > 0. 

Therefore for B^ = {i,j){i',j') — {i,j'){i',j) G Bo{S), we have Y — B^ G and we note 



that \\X — Y + B 111 < ||X — y||i regardless of the signs of Xij — yij and Xj/j/ 



If 



Xij < yij and Xi'j' < yi'j', 



IX -r 



On the other hand, if Xij > yij or Xi'ji 
Figure [2]-(i) or (ii), we have 



-B-\\i-- 
> Vi'j', 



= ll^-^'lli- 

i.e. X — Y has the pattern of signs as in 



\X -Y + B- 



\l<\\X-Y\\^. (4) 

In this case = 0, = 1 and B^ = B~ satisfy ([3]). By interchanging the role of X 
and Y, we can see that the patterns in (i) and (ii) in Figure [1] are interchangeable. Hence 
similar argument can be done for the patterns (ii) in Figure [1] and (iii), (iv) in Figure [2l 

Denote Z = Zq = X — Y . For a sequence of basic moves Bf G Bo{S), t = 1, . . . , r 
denote Zt = X — Y + Bi + -- - + Bt, t = 1, . . . , r. Based on the above arguments, we 
obtain the following lemma. The proof is easy and omitted. 

Lemma 1. ||2'||i can be reduced by Bo{S) if there exist r > and a sequence of basic 
moves Bt G Bo{S), t = 1, . . . , r such that Zt, t = . . . , r — 1, have either of the sign 
patterns in Figure [I] and Z^^ has either of the patterns in Figure [H 

This lemma will be repeatedly used from Section 13.31 on. 
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3 A necessary and sufficient condition 

In this section we give a necessary and sufficient condition on the subtable sum problem 
so that a Markov basis consists of basic moves, i.e. Bq{S) forms a Markov basis for As- 
Figure [3] shows patterns of 2 x 3 and 3x2 tables. A shaded area shows a cell belonging 
to S. Henceforth let a shaded area represent a cell belonging to S or rectangular blocks 
of cell belonging to S. We call these two patterns in Figure [3] the pattern V and P*, 
respectively. Then a necessary and sufficient condition is expressed as follows. 



Figure 3: The pattern V and 



Theorem 1. Bo{S) is a Markov basis for As if and only if there exist no patterns of the 
form V or "P* in any 2x3 and 3x2 subtable of S or S'^ after any interchange of rows 
and columns. 

We give a proof of Theorem [1] in the following subsections. Note that if Bq{S) is a 
Markov basis for As, then it is the unique minimal Markov basis, because the basic moves 
in Bo{S) are all indispensable. 
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The outline of this section is as follows. Section 13.11 gives a proof of the necessary 
condition. In Section 13.21 we introduce two patterns of 5, 2 x 2 block diagonal set and 
triangular set, and show that S or S'^ contain patterns of the form V or "P* if and only if 
S is equivalent to either of the two patterns. Then the sufficiency can be rewritten that 
Bq{S) forms a Markov basis for S equivalent to a 2 x 2 block diagonal set or a triangular 
set. In Section [3.31 we prepare some ingredients to prove the sufficiency. In Section [37il 
and Section [331 we show proofs of the sufficient condition for 2x2 block diagonal set and 
triangular set, respectively. 



3.1 A proof of the necessary condition 

The necessary condition of Theorem [1] is easy to prove. 

Proposition 1. If S or contains the pattern V orV^ , Bo{S) is not a Markov basis for 
As- 

Proof. Assume that 5* has the pattern V. Without loss of generality we can assume that 
V belongs to {{i,j) | z = 1, 2, j = 1, 2, 3}. Consider a fiber such that 

• xi+ = X2+ = 2, x+1 = x+2 = 1, x+3 = 2; 

• Xi+ = and x+j = for all (z,j) ^ {{i,j) M = 1,2, j = 1,2,3}; 



E 



Then it is easy to check that this fiber has only two elements shown in Figure HI Hence 
the difference of these two tables 



B 



1 


1 


-2 


-1 


-1 


2 



(5) 



is an indispensable move. Therefore if S has the pattern V, there does not exist a Markov 
basis consisting of basic moves. When S has the pattern P*, a proof is similar. 

□ 



1 1 



1 1 



Figure 4: Two elements of the fiber 



It is of interest to note that the toric ideal for the 2x3 table with the pattern V 
of S* is a principal ideal generated by a single binomial corresponding to ([5]) whose both 
monomials are non-square-free. 
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3.2 Block diagonal sets and triangular sets 

After an appropriate interchange of rows and columns, if S satisfies that 

•S* = \ i<r,j <c}U {(z, j) \ i> r,j > c} 

for some r < R and c < C, we say that S is equivalent to a 2 x 2 block diagonal set. 
Figure [5] shows a 2 x 2 block diagonal set. A 2 x 2 block diagonal set is decomposed into 
four blocks consisting of one or more cells. We index each of the four blocks as in Figure 
O Note that S" is a 2 x 2 block diagonal set if and only if S''^ is a 2 x 2 block diagonal set. 

For a row index i, let J^{i) = G S} denote a slice of S at row i. If for every 

pair i and i', either J'{i) is a subset of J^{i') or J'{i) is a superset of we say that 

S is equivalent to a triangular set. A triangular set is expressed as in Figure E] after an 
appropriate interchange of rows and columns. In general, if we allow transposition of 
tables, triangular sets can be decomposed into nx(n + l)ornxn blocks as in Figure 
[6l Figure [6] shows examples of n x (n + 1) and n x n triangular sets with n = 4. We 
index each block as in Figure [61 Let JF^ be a fiber of an x [n + 1) triangular set. 
Define J'n+i = {j \ € Ik,n+i}- Then we note that if JF^ satisfies J2f=i^ij = 
j £ J'n+i, the fiber is equivalent to a fiber for an n x n triangular set. Hence an n x n 
triangular set is interpreted as a special case of an n x (n + 1) triangular set. Hereafter 
we consider only n x (n + 1) triangular sets and let a triangular set mean n x (n + 1) 
triangular set. Note also that S* is a triangular set if and only if S*^ is a triangular set. In 
other words, a triangular set is symmetric with respect to 180° rotation of the table. 

Proposition 2. There exist no patterns of the form V or in any 2x3 and 3x2 

subtable of S after any interchange of rows and columns if and only if S is equivalent to 
a 2 X 2 block diagonal set or a triangular set. 

Proof. Assume that S does not contain V and and that S contains a (cell-wise) 2x2 
crossing sub-pattern presented in Figure [3 Without loss of generality the crossing pattern 
belongs to {(i, j) | i = 1, 2, j = 1, 2}. Since S does not contain V and P*, {(i, j) | i = 1, 2} 
and I j = 1,2} have to have the pattern as in Figure [8] after an appropriate 



: : : : :Ji:i: : : 


2:12 


X21 





Figure 5: 2 x 2 block diagonal set 
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'.Tri '.'i'. '. T.t'a'. 



Il5 



'.Tri'X.'.TM- 



l2i 1 125 



.l2i;.i.;.;X>. 



L33 I X34 



■'Zs.i.'-.'.Zzi- 



■ :Zu: ■ 



1i2 I 143 1 Iaa 



1i2 I 243 



Figure 6: (block-wise) 4x5 and 4x4 triangular sets 



Figure 7: 2x2 (cell-wise) crossing pattern 



Figure 8: The pattern of {(i, j) | i = 1,2} and {(i, j) | j = 1,2} 



Figure 9: The pattern of S which has a (cell-wise) 2x2 crossing pattern 
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interchange of rows and columns. In the same way the rest of the table | ^ > 3, j > 

3} has to have the pattern as in Figure [H It is clear that the pattern in Figure [9] is 
equivalent to a 2 x 2 block diagonal pattern after interchanging rows and columns. 

From the definition of triangular set, S is not equivalent to a triangular set if and only 
if there exists i, i', i ^ i', and j, j', j ^ j', such that j G J7'(z), j ^ J{i')^ f ^ and 
i' ^ J{i)- But this is equivalent to the existence of a 2 x 2 crossing pattern. □ 

3.3 Signs of blocks 

Based on Proposition [2], for the sufficient condition of Theorem [1] we only need to show 
that Bo{S) forms a Markov basis for S equivalent to a 2 x 2 block diagonal set or a 
triangular set. As mentioned above, a 2 x 2 block diagonal set and a triangular set can 
be decomposed into some rectangular blocks. In general each block consists of more than 
one cell. For the rest of this section, we use the following lemma. 

Lemma 2. Assume that S is equivalent to a 2 x 2 block diagonal set or a triangular 
set. Suppose that Z = {zij}(^ij)^x contains a block Im such that {i,j) G Iki, {i'jj') G ^ki, 
Zij > and Zi/ji < 0. Then \\Z\\i can be reduced by Bq{S). 

Proof. Suppose that j = j' and i ^ i' . Since any row sum of Z is zero, there exists j" 
such that Ziyi < as presented in Figure [TOl Hence Z contains the sign pattern of Figure 
[T]-(i) and \Z\x can be reduced by i3o(S'). When i = i' and j ^ j', we can show that 
can be reduced by Bi^{S^ in the similar way. 

Next we consider the case where i ^ i' and j ^ j'. If Zi'j 7^ or Zijr 7^ 0, we can 
reduce ||^||i by using the above argument regardless of the signs of them. So we suppose 
Zi'j = and Zij' = 0. There exists j" such that Zij" < as presented in Figure [TTI-(i). 
If {i,j") G 2kh we can reduce \\Z\\i by using the above argument. If {i,j") ^ X^/, let 
B = {i,j"){i',f) - {ij'){i',f) G Bo{S) and let Z' = {4} = Z + B. Since Zij. < and 
Zi'j' < 0, we have \\Z'\\i < \\Z\\i. We also have z^j > and < 0. Since (i, j), {i,j') G Im, 
\\Z'\\i can be reduced by Bo{S). Therefore Z satisfies the condition of Lemma[T]and 
can be reduced by Bo{S). □ 





* 







Figure 10: Z when j = j' 
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fiii) Z' 



Figure 11: Z and Z'] when j 7^ j' 
'0— " represents the cell which is nonpositive) 



Let "0+" and "0—" represent the cell which is nonnegative and nonpositive, respec- 
tively, as in Figure [TTJ 

From Lemma [2] in order to prove the sufficient condition of Theorem [H we only need 
to consider the case where Z does not have a block with both positive and negative cells. 
If all cells in Xu are zeros, we denote it by = 0. If Xm 7^ and all nonzero cells in Xki 
are positive, we denote it by Xm > and we say Xki is positive. We define Xm < 0, Xm > 
and Xki < in the similar way. Then we say X^i is negative, nonnegative and nonpositive, 
respectively. Then we obtain the following lemma. 

Lemma 3. Assume that S is equivalent to a triangular set. Suppose that Z has four 
blocks Xki, Ik'i, Iki' (in-d Xk>i> which have either of the patterns of signs as follows, 

II' I V II' II' 

k + — k * — k — + k * + 

k' - * k' - + k' + * k' + - 

(i) (ii) (iii) (iv) 

where * represents that the sign of the block is arbitrary. If there exist i,i',j,j' such that 
e Xki, (i'J) e Xk'i, (ij') e Xku {i'J') e Xk'u, and B{i,i';jJ') = {iJ){i'J') - 
{i,j'){i',j) G Bo{S), then \\Z\\i can be reduced by Bo{S). 

Proof. Assume that the four blocks have the pattern of signs (i) and that 



< 0, (ibjb) e Xk'i, Zi^j^ < 0, {ic,jc) e ^ki'- 



We note that {icjb) e Iki- If Zi^^ > OT Zi^j^ > 0, ll-Z^lli can be reduced by i3o(>S'). Suppose 
0. Let B = (ife, jc) - iibJc){icJb)- Denote Z' = {4}(ij)ei = Z + B. 



Then we have 2^'^^-^ < 0. Since there exists {iaiJa) ^ ^ki such that z'^^^^^ > 0, Z' has both 
positive and negative cells in Xki- Hence can be reduced by BqIS) from Lemma [H 

Proofs for the other patterns are the same by symmetry. □ 
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Ja Jh Jc Ja Jh Jc Ja Jb Jc 

















0- 



































(i) Z (ii) B (iii) Z' 

Figure 12: Z and Z' when Zi^j^ = Zi^j^ = 

3.4 The sufficient condition for 2x2 block diagonal sets 

In this subsection we give a proof of the sufficient condition of Theorem [T] when S is 
equivalent to a 2 x 2 block diagonal set. 

Proposition 3. If S is equivalent to a 2 x 2 block diagonal set, Bo{S) is a Markov basis 
for As. 

Proof. Suppose that Z ^ 0. li Z contains a block Xki which has both positive and negative 
cells, ||2'||i can be reduced by Bq{S) from Lemma [H 

Next we suppose that all four blocks are nonnegative or nonpositive. Without loss of 
generality we can assume that Xn > 0. Since all row sums and column sums of Z are 
zeros, we have X12 < 0, X21 < and X22 > 0. On the other hand, since X](ij)e5% ~ ^' 
we have X22 < 0. However this implies Z = and contradicts the assumption. □ 



3.5 The sufficient condition for triangular sets 

In this subsection we give a proof of the sufficient condition of Theorem [1] when S is 
equivalent to an n x (n + 1) triangular set in Figure [61 We only need to consider this case 
if we allow transposition of the tables and because of the fact that an n x n triangular set 
can be considered as a special case of an x (n + l) triangular set as discussed in Section 

In general, as mentioned, each block consists of more than one cell. However for 
simplicity we first consider the case where every block consists of one cell. As seen at the 
end of this section, actually it is easy to prove the sufficient condition of Theorem [1] for 
general triangular set, once it is proved for the triangular sets with each block consisting 
of one cell. Therefore the main result of this section is the following proposition. 

Proposition 4. Suppose that S is equivalent to ann x {n + 1) triangular set in FigurelE 
and every block consists of one cell. Then B(){S) is a Markov basis for As. 

We prove this proposition based on a series of lemmas. In all lemmas we assume that 
S is equivalent to an n x [n + 1) triangular set. If n = 1, each fiber has only one element. 
Hence we assume that n > 2. 



12 



Lemma 4. // Z contains a row ia such that the signs of Zi^i and Zi^^n+i a-fc different, 
then \\Z\\i can he reduced by Bq{S). 

Proof. Without loss of generality we can assume that Zi^i > and Zi^^^+i < 0. Since 
Y17=i ^i,n+i = 0, there exists ib such that Zj^^n+i > as presented in Figure [131 Hence if 
we set B = {ia,n + l){ib,l) - {ia,l){ib,ri + 1), B e Bo{S) and \\Z + B\\i < \\B\\i. □ 




Figure 13: The case of n = 3, = 1 and = 3 



Lemma 5. Suppose that Z has three rows ia < ib < ic satisfying either of the following 
conditions, 

(i) Zi^i > 0, Zi^i < and Zi^i > 0; 

(ii) Zi^i < 0, Zi^i > and Zi^i < 0; 
Then \\Z\\i can he reduced hy Bo{S). 

Proof. It suffices to prove the case of (i). Since Zi^i < 0, there exists j such that 2 < 
j < n + 1 and Zi^j > 0. If {ib,j) E S as presented in Figure [14]- (i), B = {ia, j){ibA) ~ 
{ia,l){ib, j) G Bo{S) and + < ||2'||i. If {ib,j) ^ 5 as presented in Figure [14]- (ii), 
B' = {ic,j){ib, 1) - {tc, l)(^6, j) G BoiS) and \\Z + B'\U < \\Z\\,. □ 



lb 



+ 



(i) {^b.3)^S (ii) {ib,3)iS 

Figure 14: The case of n = 3, ia = 1, 4 = 2 and = 3 
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Lemma 6. Suppose that Z contains four rows ia, %, ic and id satisfying 

Zi^i > 0, Zi^^n+i > 0, Zi^i < and Zi^^n+i < 
and satisfying either of the following conditions, 

(i) ia<ic< ib, (i') ia<id< (ii) ib < ic < ia, (ii') ib < id < ia ■ 
Then \\Z\\i can be reduced by Bo{S). 

Proof. Suppose (i) ia < ic < ib- Since any row sum is zero, there exists j such that 
Zi^j > and j > 2. If (zc, j) eS oi j = n + 1, B = {ia,j){ic, 1) - {ia, l)(«c, j) e Bo{S) and 
||Z + < 11^ 111 (Figure [T5] shows an example for this case). Suppose that {ic,j) G 
and j ^ n + 1. If 2i,,n+i < 0, B = {ib,j){ic,n + 1) - {ib,n + l){ic,j) e ;Bo('S') and 
||Z + i?||i < ||Z||i. If = 0, llZ'lli = ||Z + i?||i < ||Z||i and z^^^, > 0. Since 

2;-^ 1 > and z'^^ ^ < 0, ||Zi|| can be reduced by Bo{S). Hence ||^||i can be also reduced 
by Bo{S) (Figure [T6] shows an example for this case). 



+ 



(i)(^c,j)eS (ii)j=n+l 
Figure 15: The case of n = 3 and ia < ic < ib 



+ 



-1 



0+ 



(i) Z (ii) B (iii) Z' 

Figure 16: The case of n = 3 and ia < ic < ib 

In the case (i') ia < id < ib, we can prove the lemma in the same way by the symmetry 
of n X [n + 1) triangular pattern. 

Suppose that (ii) ib < ic < ia or (ii') ib < id < ia- If ^i^.n+i < or Zi^^i < 0, 
the lemma holds from Lemma HI So we suppose that Zi^^n+i > and Zi^^i > 0. Let 
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-1 



+ 



(i) Z (ii) B (iii) Z' 

Figure 17: The case of n = 3 and ib < ic < ia 



B = [ia, n + l){ibi 1) — {iai n + 1). Then we have ||Z'|| = ||Z + < ||^||i and 
^'ia,n+i > 0' > 0- Since z'^^ i = Zi^^i < and z[^^^_^_i = Zi^^n+i < 0, we can prove the 
lemma by applying the above argument (Figure fTTl shows an example for this case). □ 

From the definition of Li-norm, ||^||i can be reduced by ^30(5"), if and only if || — Z\\i 
can be reduced by Bo{S). Thus without loss of generality we can assume that 2:11 > 0. 
From Lemma m [S] and El it suffices to show that can be reduced by Bo{S) if Z 

satisfies 

3zo > 1 s.t. 

Zii > and Zi^n+i > for i < io, (6) 
Zii < and 2i,n+i < for i > io, 

as shown in Figure | 

















;£):^;:;:;:;:;:;:;: 


i 0- 
4 1 ; 






i i : 0- 

1 1 ' 
1 1 



Figure 18: The case of n = 4 and ^o = 2 



We look at rows of Figure [18] from the bottom and find the last row ii such that at 
least one of Zi^^i or is negative, i.e., define ii by the following conditions. 

(i) Zi^i<0 and Zi^^n+i < 0; 

(n) Zi^i < or 2ii,n+i < 0; (7) 
(ni) Zii = and Zi^n+i = for i > ii; 

Note that if there exists no ii satisfying these conditions, then the first column and the last 
column of the table consists of only zeros and we can use the induction on n. Therefore 
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for Lemmas [3-[9] below, we assume that there exists ii satisfying ([7j). We also note that 
ii > if) when ii exists. 



Lemma 7. Suppose Z satisfies ([^ and define ii by ^ assuming that ii exists. \\Z\\i can 
he reduced if there exists Zy > for some {i,j) G 5* and i > ii- 

Proof. Consider the case Zi-^i < 0. We note that there has to exist i' < iq such that Ziii > 
from the condition Q). Suppose that i = ii. Let B he B = (ii,l)(z',j) — (zi, 1). 
Then \\Z + B\\i < \\Z\\i. Suppose that i > ii and Zi^j < for {ii,j) G S. There has to 
exist / such that Zi-^ji > 0. Let B = (zi, j)(z, j') — (zi, j). If Zi^^j < or Ziji < 0, 
llZ'lli = \\Z + B\\i < \\Z\\i. If Zi,j = or Zij- = 0, \\Z'\\i = \\Z + 5||i < \\Z\\i. As shown 
in Figure fT9| since z[,i > 0, z'^^-^ < and z'^^j > 0, \\Z'\\ can be reduced by Bo{S). Hence 
II Z 111 can also be reduced by Bo{S). 

Next we consider the case Zi-^^n+i < 0. Then there has to exist i' < such that 
Zii^n+i > from the condition ([6]). When i = ii, let B = {ii,n + l){i',j) — (ii, + 1). 

Then \\Z + i?||i < ||Z||i. When i > ii and Zi-^j < for {ii,j) G S", a similar proof to the 
case < can be given as shown in Figure [201 

□ 
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Figure 19: The case of n = 6, {i(),ii,i,i' , j, j') = (3,4,5,2,2,4) and Zi-^^i < 
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Figure 20: The case of n = 6, {io,ii,i,i' , j, j') = (3,4,5,2,2,4) and Zi^^n+i < 
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We define some more sets. Let S"^ and S^, i = 2, . . . ,n, he the sub-triangular set of 
defined as 

= {it',f) eS^\f^n + 1}, = {(z',/) eS^\t'<t, f ^ n + 1}, 
respectively. Figure fill shows S'^ and S'f for n = A, i = 4. We note that 

for all Z, because the last column sum is zero and J2{ij)es'' % ~ ^- define S^, 

z = 2, . . . , n, by 




Figure 21: S*^ and for = 4 and z = 4 

From Lemma [7] it suffices to consider Z such that Zij < for all G S and i > ii- 
The following lemma states a property of such a Z. 

Lemma 8. Suppose that Z satisfies ^ and define ii by ^ assuming that ii exists. 
Furthermore assume that Zij < for all {i,j) & S , i > ii. Then 

Proof. Assume that 

{^.i)e5- 

Since the roles of Zi-^i and -Zii,n+i are interchangeable, we assume that \zi-^^i\ > 0. Then 
there exist nonnegative integers w},, wZ^^ and the set of cells S' C S~ and S" C 5," 
satisfying 

(ij)e5' (ij)6S" (ij)es- 

17 



< I 1 1 ' < ,n+i I • 

iLj)eS' {i,j)es" 

S' and S" may have overlap if \zij\ > 2 for some cell (z,j) G S'~. For (z, j) G 5^^, let 
B^{i,j) and B^~^^{i,j) be defined by 

j) = (z, 1) - l)(^i, j), j) = J)(^i,^ + 1) -{^,n + l)(zi, j)- 

We note that B\i,j) G i3o(5) and B''+^{i,j) G So(5) for any (z,j) G Denote 

^' = {4}(M)ex = ^+ 5^ wl^B\z,j)+ Y <;'B-^\z,j). (10) 

(i,i)65' (i,j)&S" 

Then we have z^'^^ < 0, < 0, and z[j > for all (z,j) G S'j'^^. This implies that 

E 4<o- (11) 

On the other hand from the condition of Lemma [8] 

n n+1 n 

E 4 = EE4-( E 4 + E<"+i)>0' 

which contradicts (fTT!) (See Figure [22|) . □ 
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Figure 22: The case of = 6 and (zq, zi) = (3, 5) 

Lemma 9. Suppose that Z satisfies ^ and the conditions of Lemmal^ Then 

(i) ii > 3; 

(ii) If ii = 3, \\Z\\i can be reduced. 
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Proof, (i) It is obvious that ii > 2. Suppose ii = 2. Since any row sum of Z is zero, we 
have 

from (ii) and (iii) of ([7j). However this contradicts ([8j). 

(ii) When ii = 3, S3 = {(2, n)} and Z2n < from Lemma [H If Z21 > 0, we have 2:31 < 
from (iii) of ([7]). Therefore S = (3, 1)(2, n) - (2, 1)(3, n) satisfies + < ||Z||i (Figure 
[23]-(i)). If 22,n+i > 0, ^3,„+i < Ofrom (iii) of dZD. Hence 5 = (2, n)(3, n+l)-(2, n+l)(3, n) 
satisfies \\Z + i?||i < \\Z\\-\ (Figure [23]- (ii)). Next we consider the case of Z21 = Z2,n+i = 0. 
Then zu > or Zi,„+i > 0. Suppose that zu > 0. This implies z^i < 0. Since Z2n < 0, 
there exists 2 < j < n — 1 such that Z2j > 0. Then if we set B = (1, j)(2, 1) — (1, 1)(2, j), 



we have 11 Z' 



\Z + B\U < \\Z 



|i and z'21 > 0, 



^31 



< and 4„ < (Figure [M]). Then 



||Z'||i can be reduced by l3o{S). Therefore \\Z\\i can be also reduced by Bo{S). When 
zi^n+i > 0, a proof is similar. □ 
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::Q+:: 
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(i) 2:21 > (ii) Z25 > 

Figure 23: The case of n = 4 and Z21 > or 2:25 > 
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(i) Z (ii) B (iii) Z' 

Figure 24: The case of n = 4, 2:21 = and 225 = 
By using Lemmas HHHl we give a proof of Proposition m 
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Proof of Proposition We prove this proposition by the induction on the number of rows 
n. Suppose that n = 2. Then 

zij + Z2j = for j = 1, 2, 3. 

Zll + Zi2 + Z2l = 0, Z22 + Zi3 + Z23 = 0. 

Hence 2:12 = Z22 = 0. Therefore Z is equivalent to a move in the 2x2 pattern as in Figure 
[25] with fixed row sums and column sums. It is easy to see that that proposition holds for 
this pattern. 



Figure 25: A 2 x 2 pattern 

Suppose n > 2 and assume that this proposition holds for triangular sets smaller than 
n X (n + 1). From the results of Lemmas HHHl it suffices to show that if Z satisfies (EI) 
and the conditions of Lemma [HI \\Z\\i can be reduced by Bo{S). We prove this by the 
induction on ii. 

Suppose that il > 3 and assume that Z with ii < il can be reduced by Bo{S). From 
Lemma [HI ([9]) holds. Thus there exist nonnegative integers wjj, w^^'^ and the set of cells 
S' C S~ and S" C S^^ satisfying 

wlj + <\zij\, = I 1 1 ' = I '"+1 1 • 

Let Z' be defined as in f [TO|) . Then we have \\Z'\\i < \\Z\\i. If \\Z'\\i < ||^||i, this 
proposition holds. Suppose that \\Z'\\i = \\Z\\i. Then Z' satisfies either of the following 
three conditions, 

(i) z'^^ = and 2- „+i = for i = 1, . . . , n; 

(ii) there exists i such that z'--^ 7^ or ^' does not satisfy ([6]). 

(iii) there exists i such that z[^ 7^ or 7^ ^^"^ ^' satisfies ([6]). 

In the case of (i), can be reduced by Bq{S^ from the inductive assumption on n. 

In the case of (ii), can be reduced by by Lemma [B In the case of (iii), noting that 

z'^^y = and ^^'n+i = for z > zi, \\Z'\\i can be reduced from the inductive assumption on 
il. ' □ 

So far we have given a proof when every block has only one cell. It remains to extend 
Proposition [H to general triangular sets. Based on the results of Lemma [1] and [2], we see 
that Proposition [H can be extended to the case where n>2. Then it suffices to consider 
the case of n = 1 as in Figure [23 
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Figure 26: 1x2 triangular pattern 

Lemma 10. Suppose that S is equivalent to a 1 x 2 triangular set. Then Bq{S) is a 
Markov basis for As- 

Proof. Since Z satisfies J2{ij)es^ii ~ ^ j)es'' ^ij = 0, Z has to contain both 

positive and negative cells in at least one of Xn and X12. Hence ||2'||i can be reduced from 
Lemma [H □ 

Now we have completed a proof of the sufficient condition of Theorem [1] for general 
triangular set 5*. 



4 Concluding remarks 

In this paper we consider Markov bases consisting of square-free moves of degree two for 
two-way subtable sum problems. We gave a necessary and sufficient condition for the 
existence of a Markov basis consisting of square- free moves of degree two. 

From our results, if S contains a pattern "P or P*, a Markov basis has to include a move 
with degree higher than or equal to four. From theoretical viewpoint, it is interesting to 
study the structure of Markov bases for such cases. Our results may give insights into the 
problem. However it seems difficult at this point and left to our future research. 

Consider a particular fiber with x{S) = in the subtable sum problems. Then Xij = 
for all {i,j) G S. This implies that this fiber is also a fiber for a problem where all cells 
of 5" are structural zeros. Therefore Markov bases for the subtable sum problems for S 
are also Markov bases for a problem where all cells of 5* are structural zeros. Various 
properties of Markov bases are known for structural zero problems. It is of interest to 
investigate which properties of Markov bases for structural zero problem for S can be 
generalized to subtable sum problem for 5*. 

Ohsugi ar id Hibi have been inves t igating properties of Grobner bases arising fro m 
finite graphs flOhsugi and Hibil |l999al | , lOhsugi and Hibil |1999bl | , lOhsugi and Hibi j2005| l 



With bipartite graphs, their problem is equivalent to two-way contingency tables with 
structural zeros. From the viewpoint of graphs of Ohsugi and Hibi, the subtable sum 
problem corresponds to a complete bipartite graph with two kinks of edges. It would be 
also very interesting to investigate subtable sum problem from the viewpoint of Grobner 
bases. 

We used the norm reduction argument to prove that Bo{S) is a Markov basis. It should 
be noted that Bq{S) for the subtable sum problem is not necessarily 1-norm reducing in 
one step, even when Bo{S) is the unique minimal Markov basis. Therefore the subtable 
sum problem is worth to be considered from the viewpoint of norm reduction by a Markov 
basis. 
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